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We discuss recent results on the three-loop (next-to-next-to-leading order, NNLO) time-like splitting functions 
of QCD and the two-loop (NNLO) coefficient functions in one-particle inclusive e^e~-annihilation. These results 
form the basis for extracting fragmentation functions for light and heavy flavors with NNLO accuracy that will 
be needed at the LHC and ILC. The two-loop calculations have been performed in Mellin space based on a new 
method, the main features of which we also describe briefly. 



1. Introduction 

Observables with identified hadron(s) in the fi- 
nal state depend on the fragmentation function(s) 
for these hadron(s). A prominent example for the 
physics at hadron colliders is the p^-spectrum of 
6-flavored mesons [1-3] , 
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Unlike the parton distributions (pdf's) fi, the 
scale dependence (evolution) of the fragmentation 
functions -Dfc/B is in terms of time-like splitting 
functions. As is well-known [4-6], starting from 
two loops (i.e., NLO) the time- like and the space- 
like splitting functions of QCD differ. 

To derive the time-like splitting functions to 
any order in perturbation theory, one can ex- 
ploit the universal factorization property of the 
coUinear singularities in some particular process, 
for instance in e''"e~-annihilation to hadrons. 



f{p)+X., 



(2) 



where /(p) can be any massless on-shell parton. 
The differential cross-sections for the reaction (2) 
are known to two loops [7-10], permitting the ex- 
traction of the time-like NLO splitting functions 
and NNLO coefficient functions. 
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The full three-loop calculation, needed for ob- 
taining the NNLO splitting functions is a problem 
of remarkable complexity due to, among other 
reasons, the incomplete inclusiveness of the fi- 
nal state in this process. In contrast to the re- 
cent calculation of deep-inelastic scattering (DIS) 
[11-14], one cannot apply the optical theorem 
here. Progress in this direction is being made [15] 
which will be discussed in Sec. 3. Before we 
present in Sec. 2 an alternative approach and first 
results for the time-like NNLO splitting functions. 



2. Time-like NNLO splitting functions 

The essence of our approach is to devise a re- 
lation between the space- and time-like splitting 
functions in QCD capable of predicting the quan- 
tities with time-like kinematics from the known 
space-like ones [11,12]. This approach is a sub- 
ject with interesting history starting with the 
Gribov-Lipatov relation at leading order [16]. 
The NLO calculations [4-6] explicitly demon- 
strated the breaking of the Gribov-Lipatov rela- 
tion beyond the leading order, see also Ref. [17]. 
The interest in the subject has been re-emerging 
throughout the years [18-20], and the space- 
or time-like two-loop results were independently 
confirmed by several groups [10,21-23]. New re- 
sults beyond two loops, however, were obtained 
only very recently [24] . 

Our starting point is the factorization of mass- 
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less differential cross-sections. Specifically, we 
consider in parallel the following two reactions: 



DIS : f{p)+r ^ X, 



(3) 



and the corresponding parton level observables 

X = 

dx 



da —q" 9 



da 

dz 
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(4) 



The above differential distributions contain sin- 
gularities of coUinear origin. These, however, 
completely factorize and can be absorbed into the 
definition of the corresponding non-perturbative 
parton distribution or fragmentation functions. 
Thus, they are completely specified in terms of 
the kernels of the evolution equation. 

The observables in Eq. (3) (including Z-boson 
exchange) are expressed in terms of three struc- 
ture functions. In DIS these are known as Fi , 
and while their counterparts in e^e~ are de- 
noted as Ft, Fl and Fa- In the following we will 
consider only the relationship between the func- 
tions Fi and Ft, the others being similar. 

We first take the bare (i.e. collinear singulari- 
ties are not factored out) structure function F^' 
expanded in terms of the bare strong coupling a^, 



F^ia^^Q-") = 
6{l-x) + 



(5) 



The terms in the above series are written as 

F^\ =. 2TiS{l-x)+ni (6) 
Fi^;^ = 2J-2(5(1-.t) + (.Fi)'(5(1-x) 

+2J■l7^l +7^2 

F^^^s = 2T3S{l-x)+2TiT2S{l-x) 

+{2T2 + (Tif)ni + 2Tin2 + na . 

This iterative decomposition into form-factor (J^) 
and real-emission (7?.) parts is done according to 
the e-dependence in the soft limit, see Ref. [25]. It 



represent the closest approach to a full decompo- 
sition into contributions from the separate phys- 
ical cuts possible on the basis of Refs. [11-14]. 

To predict the time-like quantities from their 
space-like counterparts, we subject Eq. (6) to the 
following analytical continuation: first, as can be 
easily seen from the definitions (4) of the kinemat- 
ical variables in the two processes, one changes 



X ^ 1/x and 



-q . Accounting for the dif- 



ference in the two phase-space measures, one has 
to multiply the time-like expression by an overall 
factor of z^"^*^ (see [10] for more details). The 
only subtle point in the analytic continuations 
is the treatment of logarithmic singularities for 
X ^ 1, cf. Ref. [17], starting with 



ln(l — a;) ln(l — 2) — In^; -|- Itt . 



(7) 



After the continuation has been performed, we 
keep only the real parts of the continued functions 
TZ I and then re-assemble the resulting expression 
analogously to Eq. (6). This later result is iden- 
tified with the corresponding structure function 
F^ in e+e" which, according to the mass factor- 



ization, takes the form 
F^, = -e-ipW + 



-h ea, 



(1) 



(8) 



+ e26« + e^^) + ... 



and correspondingly for the coefficients Fj, ^ and 
(see Ref. [24] for their explicit expressions). 
The same procedure applies to the other two 
structure functions F^ and F\. 

The above described analytic continuation can 
in principle be applied to any order in the strong 
coupling and to any order in the expansion in e. 
Indeed, by comparing to the known expressions at 
two loops [4-10] for the time-like splitting func- 
tions P^^' as well as the coefficient functions (P"^ 
and the new a*'^-' terms, wc obscirvc; complete 
agreement - which, beyond order e^^, depends 
on the inclusion of the Itt term in Eq. (7). More- 
over, all terms proportional to e*^. k < —2 at three 
loops are also correctly predicted by the analytic 
continuation carried out in this manner. 

Despite the highly non-trivial predictions men- 
tioned above, we would like to stress the fol- 
lowing point. Unlike the per-diagram treatment 
of Refs. [4-6,17], our continuation relies on the 
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space-like evaluation of Refs. [11-14] based on the 
optical theorem. Therefore, we have incomplete 
information on the separate contributions from 
the various physical cuts, recall that Eq. (6) does 
not represent a full decomposition according to 
the number of emitted partons. Thus one has to 
be prepared for some problem at the third order, 
especially in the abelian (Cj) piece, related to tt^ 
contributions originating from phase space inte- 
grations over unresolved regions. 

There are two checks on the predicted three- 
loop time-like splitting functions: their soft be- 
havior a; ^ 1 and the vanishing of a first moment 
N = 1. The latter check shows a mismatch in a 
term of the type w'^Cpil + x'^)/{l - x) In^ x, i.e., 
the analytic continuation returns an incorrect co- 
efficient for this term which has been restored by 
the N = 1 constraint. 

Clearly, one needs to devise a second and in- 
dependent confirmation of this corrected predic- 
tion. For this purpose wc adopt the approach 
of Dokshitzer, Marchesini and Salam [26]. There 
the evolution equations for either parton distribu- 
tions or fragmentation functions /J^*^ arc rewritten 
in a unified manner for space- like (a = —1) and 
time-like (cr = -|-1) kinematics. 



d 



din 02 

^ dz 
z 



fr{x,Q^) = (9) 



with the universal splitting functions i^^fv 
sumed to describe both the time-like and space- 
like cases. 

One can easily work out the perturbative ex- 
pansion of Eq. (9). The resulting expression 
for the difference P}l\'"'ix) - (x) at NLO 

agrees with Ref. [4], while the NNLO prediction 
coincides with the (N = 1 corrected) result from 
our analytical continuation. It can be cast in the 
following very compact form with (8) denoting the 
usual convolution, 

<5P(2)?(^) = 2{[ln:r-p(i)«] (10) 



for the non-singlet kernels (^ = +, — , v) with 

2P(")€(a.) = P}Z\Hx)+P}Zl{{x) . (11) 

Obviously Eq. (10) predicts the correct vanish- 
ing first moment for P^Ji {x) and, moreover, we 
can even make a prediction for the fourth-order 
(N^LO) difference SP of the (both unknown) 
time-like and space-like non-singlet splitting func- 
tions on this basis. The difference can also be 
written in a form similar to Eq. (10), see Ref. [24] 
for explicit results. 

3. Hadron fragmentation in e'^e~ at NNLO 

Let us now turn to our method used for di- 
rect calculations of higher-order QCD corrections 
in e"'"e~-annihilation (2). In contrast to the pre- 
ceding section where we have discussed the an- 
alytic continuation in kinematic variables, here 
we describe the derivation of the two-loop cor- 
rections to the coefficient functions for light par- 
ton production. They have been first calculated 
in Ref. [7-9], later confirmed by the analytical 
continuation as outlined above and, finally, re- 
derived through an independent direct two-loop 
calculation [10]. The aim of the latter calcula- 
tion was to check the earlier results of Rijken and 
van Neerven and to derive the terms a^^) (defined 
as in Eq. (8)). The latter terms of 0{e) were a 
useful check on predictions of the analytical con- 
tinuation at two loops, and will also be needed in 
a future calculation of the three-loop corrections. 

The calculation [10] was performed directly in 
Mellin space, following the method of Ref. [15], 
which consists of two basic steps. First, one 
performs the Mellin transform before any phase- 
space or virtual integration: 



a{N) 



dzz^-^% 
dz 



[VI) 
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lnar-P(°) 



p(i)«| 



= ^ dzz"" J dPS^"' \MfS{z - f) 

= y dPS('")]M|V^"S 

where / is some function of the external and/or 
integration momenta that defines the variable z. 
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If the measure dPS^"^' contains (5-functions, they 
can be dealt with in the standard way [27]. Prom 
Eq. (12) it is clear that the effect of the Mellin 
transform is to introduce a new propagator raised 
to an abstract power. The generalization to mul- 
tiple differential distributions involving a Mellin 
transform in multiple variables is straightforward. 

Second, one applies the reduction identities 
based on integration-by-parts (IBP) to simplify 
the right-hand side of Eq. (12) and to reduce it 
to master integrals. The fact that one of the prop- 
agators is not a fixed integer but an abstract pa- 
rameter does not present any complications. This 
is because the IBP identities can be cast in the 
usual form (where the powers of all propagators 
are fixed integers) through a shift in the variable 
corresponding to the power of the propagator /. 
The resulting recurrence relations explicitly de- 
pend on N and can be solved in a standard way 
like, e.g., the algorithm of Laporta [28] as imple- 
mented in the Maple package AIR of Ref . [29] . 

The master integrals resulting from the solu- 
tion of the IBP identities implicitly depend on 
TV, too, and one can derive for them difference 
equations in N. To solve the latter we ap- 
ply techniques used previously in the context of 
DIS [11,22,25] and perform the required expan- 
sions in e with the packages Summer [30] and 
XSUMMER [31]. By solving the difference equa- 
tions one extracts the complete dependence on 
the Mellin variable A'^ without the need for ex- 
plicit evaluations of any Feynman integral. Fi- 
nally, to fully specify the solution of a difference 
equation of degree k, one has to provide k initial 
conditions. Here, the simplest option, which we 
have used in Ref. [10], is the direct evaluation of 
the master integrals for iV = 1, . . . , k. 

All but one of the master integrals in Ref. [10] 
satisfy difference equations of first order, the 
other being of second order. As pointed out in 
Ref. [15], however, not all of the integrals speci- 
fying the initial conditions are independent, since 
the IBP reductions exhibit many additional rela- 
tions after fixing N, say to A'' = 1. Typically, for 
the initial conditions one has to explicitly eval- 
uate less objects than the number of master in- 
tegrals. In the case of Ref. [10] , we encountered 
six real-real and respectively five real- virtual mas- 



ter integrals along with seven independent initial 
conditions. 

Finally, a particularly powerful feature of the 
method of Ref. [15] is the fact that the master 
integrals depend only on the process but are in- 
dependent of the particular observable. For that 
reason all integrals needed for the initial condi- 
tions could be taken over from Ref. [32]. We con- 
firmed those results and extended them to higher 
orders in e. 

4. Energy spectrum of 6-quarks in e+e~ 

The most important application of the time- 
like QCD splitting functions is the description 
of hadron fragmentation in processes sensitive 
to coUinear radiation, like in Eqs. (1) or (2). 
For instance, the process independent fragmenta- 
tion functions D^/b of 6-flavored mesons are ex- 
tracted from measurements in e^e~-annihilation. 
Along with high-quality data, a consistent NNLO 
description of light-flavor fragmentation requires 
the knowledge of the three-loop time-like split- 
ting functions and the two-loop coefficient func- 
tions for the production of massless partons. Ad- 
ditionally, within the perturbative fragmentation 
function (PFF) formalism [33] for heavy-flavor 
fragmentation one also has to supply initial con- 
ditions now known to two loops [34,35]. Flavor 
threshold crossing conditions may also have to be 
taken into account [36,37]. 

As an application of our results we can read- 
ily predict the two-loop (NNLO) energy spectrum 
of massive 6-quarks in e+e~-annihilation based 
on the PFF formalism [33]. Now this formalism 
can be applied completely and consistently at two 
loops, since all fermion- initiated components are 
presently known [34]. The only two-loop com- 
ponent of the PFF not required in e+e" at the 
two-loop fixed-order level is the gluon-initiated 
one [35]. 

The spectrum of a massive 6-quark in e~^e~- 
annihilation, up to power corrections 0{m), 
reads: 

^{z,Q,m)= (13) 
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where is the factorization scale and Q represents 
the characteristic hard scale of the reaction. 

For the case of e"'"e~-collisions one usually takes 
Q = y/s and m to be the pole mass of the heavy 
quark. Inserting the explicit expressions for the 
coefficient functions and the PFF's one can verify 
that indeed all dependence on the factorization 
scale /i completely cancels. 

A particular feature of Eq. (13) which we would 
like to point out is the characteristic dependence 
1/2: at low z that appears for a first time 
at two loops. This is entirely due to the purc- 
singlet fermion emission and similar to the well- 
known small-a; behavior of structure functions in 
DIS [11-14]. 

Detailed phenomenological studies of the en- 
ergy spectrum dab/dz of Eq. (13) and related ob- 
servables will appear in a forthcoming publication 
[38]. 

5. Conclusions 

In these proceedings we have reviewed re- 
cently derived results for the non-singlet compo- 
nents of the three- loop (NNLO) time- like splitting 
functions and the complete two-loop (NNLO) 
partonic cross-sections for one-particle inclusive 
hadro-production in e"'"e~-annihilation. 

The former results [24] have been obtained 
by utilizing an analytic continuation in the 
proper kinematic variables in connection with the 
idea [26] of a universal splitting function gov- 
erning the space- and time-like parton evolution 
equations. 

The latter results [10] provided a check on ear- 
lier calculations by Rijkcn and van Neervcn [7 9] 
and extended their results to include higher or- 
der terms in e. The derivation of these so far 
unknown terms of 0{e) was essential in verifying 
predictions of the space-like to time-like analyt- 
ical continuation and will also be required in a 
future direct three-loop evaluation of this process. 

The derivation of the two-loop coefficient func- 
tions in e+e~ was a first application of a new 
method for evaluation of differential distributions 
directly in Mellin space [15]. This application 
demonstrated the anticipated eSiciency of the 
method both, in terms of solving the IBP identi- 



ties and, in the decrease of the number and the 
complexity of the master integrals that require 
separate treatment. 

The above results, together with the available 
two-loop contributions [34] to the PFF open the 
possibility for studies of light and heavy hadro- 
production consistently at NNLO. For example, 
the extraction of 6-fragmentation functions at 
this level of accuracy has important applications 
not only at a future ILC but also at the LHC, for 
instance in measurements of the pr-spectrum of 
i3-mesons or in precision determinations of the 
top- quark mass. 
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